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RESUMEN 

Presentamos los parametros fundamentales de enrojecimiento, metalicidad, edad y distancia de los ciimulos 
abiertos poco estudiados Be 89, Ru 135 y Be 10, derivados de la fotometria CCD UBVRI de los mismos. Las 
distancias y edades de los ciimulos se obtuvieron ajustando las isocronas apropiadas a las secuencias observadas 
de los ciimulos para cinco diferentes diagramas color-magnitud. Los promedios pesados de modulos de distancia 
y de distancias heliocentricas, (( Vq-Mv), (i(kpc)) encontrados son los siguientes: (lL m 90±0™06, 2.4±0.06) para 
Be 89, (9™58±0™07, 0.81±0.03) para Ru 135 y (ir. n 16±0™06, 1.7±0.05) para Be 10, mientras que los promedios 
pesados para las edades (log(A), A(Gyr)) son (9.58 ± 0.06, 3.8 ± 0.6) para Be 89, (9.58 ± 0.06, 3.8 ± 0.7) para 
Ru 135 y (9.06 ± 0.05, 1.08 ± 0.08) para Be 10. 

ABSTRACT 

The fundamental parameters of reddening, metallicity, age, and distance are presented for the poorly studied 
open clusters Be 89, Ru 135, and Be 10, derived from their CCD UBVRI photometry. By fitting the appropriate 
isochrones to the observed sequences of the clusters in five different color-magnitude diagrams, the weighted 
averages of distance moduli and heliocentric distances ((Vq-Mv ), ci(kpc)) are (11™90 ± 0™06, 2.4 ± 0.06) for 
Be 89, (9 m 58 ± (F07,0.81 ± 0.03) for Ru 135, and (ll m 16 ± 0™06, 1.7 ± 0.05) for Be 10, and the weighted 
averages of the ages (log(A), A(Gyr)) are (9.58 ± 0.06, 3.8 ± 0.6) for Be 89, (9.58 ± 0.06, 3.8 ± 0.7) for Ru 135, 
and (9.06 ± 0.05, 1.08 ± 0.08) for Be 10. 



1. INTRODUCTION 

Galactic open clusters, which contain a few tens 
to a few tens of thousands of stars and are a few 
parsecs across, are sparsely populated, loosely con- 
centrated, and gravitationally bound systems. With 
systematic image searches and follow-up photomet- 
ric surveys, new open clusters are currently be- 
ing discovered. By fitting the photometric observa- 
tions of open clusters to synthetic photometry result- 
ing from stellar models (i.e. theoretical isochrones), 
which include the newest input physics, stellar struc- 
ture, and differing heavy-element abundances, fun- 
damental parameters such as interstellar reddening, 
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metallicity, distance modulus, and age can be pre- 
cisely and accurately determined. These param- 
eters have great importance concerning the age- 
metallicity relation and the meta l-abundance gra- 



dient in the Galactic disk (e.g. Cameron 19851 : 



ICarraro fc Chiosi Ill994l iFriel Ill995lh as well as the 



luminosity and mass f unctions of the open clusters 
( Piskunov et al. 20081 ). Open clusters are also very 
useful for testing the stellar evolutionary models, 
given that their stars were formed at the same time, 
out of the same cloud, and under similar environ- 
mental conditions. Thus, open clusters are ideal en- 
tities for the study of stellar evolution since physi- 
cal properties are tightly constrained, being mainly 
distinguished by the stellar mass, so that theoreti- 
cal models of stellar formation and evolution can be 
compared with real clusters without excessive com- 
plications. For these analyses, the fundamental pa- 
rameters such as interstellar reddening, metallicity, 
distance modulus, and age should be determined as 
precisely and accurately as possible. 

In Galactic studies, one of the more severe obser- 
vational limitations is due to the absence of photo- 
metric data for nearly half of the approximately 1500 
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open clusters known. Furthermore, there is a lack 
of homogeneity in the observations and analyses of 
the clusters studied. The catalogue of Lynga (1987), 
that resulted from a collection of data from many 
different sources and which includes 422 open clus- 
ters, constituted the observational basis for a large 
number of astronomical studies, led to important 
conclusions about the Galactic disk, and has been 
very useful for planning subsequent observations by 
other astronomers. However, this catalogue has been 
built from parameters obtained by various authors, 
with diverse observing techniques, distinct calibra- 
tions, and different criteria for determining the stel- 
lar ages, rendering it very inhomogeneous and lim- 
ited for studies requiring precision in the measure- 
ment of these fundamental parameters. As an ex- 
ample of the precision and accuracy that one can 
expect due to the effects of these inhomogeneities, 
we refer to Janes & Adler (1982), who found that 
distance moduli of a given cluster obtained by two 
or more authors have a mean difference of 0™55. 

Within the Sierra San Pedro Martir, National As- 
tronomical Obs ervatory (hereafter SPM) open clus 



ter p roject (cf. ISchuster et al. 2007 ; Michel et al. 



2010), the aims are the following: 



1. A common UBVRI photometric scale for open 
clusters, 

2. An atlas of color-color and color-magnitude di- 
agrams for these clusters, 

3. A homogeneous set of cluster reddenings, dis- 
tances, ages, and, if possible, metallicities, 

4. An increased number of old, significantly red- 
dened, or distant, open clusters, and 

5. A selection of interesting clusters for further 
study. 

The open clusters for the present study were se- 
lected from the large (and most complete) catalogue 
, "Optically visib le open Clusters and Candidates" 
( Dias et al. 2002), which is now also available at the 
CDS El (Centre de Donnees Astronomiques de Stras- 
bourg). This work aims to provide the fundamental 
parameters of reddening, metallicity, distance mod- 
ulus and age for the open clusters Be 89, Ru 135, 
and Be 10. Our final intention is to publish a set of 
homogeneous pho tometric UBVRI data for over 300 



Gala ctic clusters (jSchuster et al. 
2010h 



2007tlTapia et al. 
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This paper is organized as follows: §2 describes 
the observations and reduction techniques. §3 con- 
tains the derivation from the UBVRI photometry of 
reddening and metallicity of the clusters from two- 
color diagrams, and the inference of distance moduli 
and ages from color-magnitude diagrams. Their un- 
certainties are also discussed. Comparisons of these 
parameters with previous results from the literature 
are made in §4, and the conclusions are given in §5. 

2. OBSERVATION AND REDUCTION 
TECHNIQUES 

2.1. The observations 

This CCD UBVRI project of northern open clus- 
ters has been undertaken at SPM using always the 
same instrumental setup (telescope, CCD-detector, 
and filters), observing procedures, reduction meth- 
ods, and system of standard stars (Landolt 1983, 
1992). A par focal set of UBVRI Johnson-Cousins 
filters was used for our observations. The 0.84- 
m f/13 Cassegrain telescope hosted the filter-wheel 
"Mexman" provided with the SITE#1 (SI003) CCD 
camera, which has a 1024x1024 square pixel array 
and a 24/imx24/im pixel size; this CCD has non- 
linearities less than 0.45% over a wide dynamical 
range, no evidence for fringing even in the I band, 
and Metachrome II and VISAR coverings to increase 
sensitivity in the blue and near ultraviolet. The sky- 
projected pixel size was 0'.'393, and the field of view 
of the detector was 6.73 x 6.73 arcmin 2 . Here the re- 
sults of UBVRI images for the open clusters Be 89, 
Ru 135, and Be 10 are presented, which were ac- 
quired in July 2001 (Be 89 and Ru 135) and Febru- 
ary 2002 (Be 10). The exposure times were typi- 
cally 3 x 240 seconds for the U filter, 3 x 180 for 
B, 3 x 100 for V, 3 x 100 for R, and 3 x 120 for 
/. Several standard-star fields from Landolt (1992) 
were observed nightly to permit the derivation of 
the photometric transformations to the Johnson- 
Cousins' system and the atmospheric extinction co- 
efficients. For the July 2001 observing run, seven 
Landolt groups were used, containing 26 different 
standard stars with color ranges, — 0^25 < (B- 
V) < +l m 14, -l m 09 < (U-B) < +l m 14, and 
-(T30 < (V-I) < +l m 14. Sixteen to twenty-five 
observations of these Landolt standards were made 
per night. For the February 2002 run, eight Lan- 
dolt groups were employed, containing 35 different 
standard stars with color ranges, —0^30 < (B- 
V) < +l m 42,-l n T8 < (U-B) < +1™27, and 
-CF28 < (V-I) < + l™77. Fifty-two to seventy-two 
observations of these Landolt standards were made 
per night, except one night cut short by clouds, when 
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Fig. 1. DSS red-filter images of the Galactic open clus- 
ters Be 89 (panel a), Ru 135 (b), and Be 10 (c). The 
regions analyzed with the ELIPSE inspection tool, to de- 
rive first estimates for the fundamental parameters, are 
enclosed by ellipses. Orientation as usual: north is up, 
and east to the left. 



only 15 observations were managed. The standard- 
star fields have been observed with exposures of 
1 x 240 seconds for the U filter, 1 x 120 for B, 1 x 60 
for V, 1 x 60 for R, and 1 x 60 for I. The observed 
Landolt fields and the number of associated stars in 
each one are summarized in Table 1. 

TABLE 1 

Landolt 's fields of standard stars 



July 2001 




February 2002 




region 


N 3tda 


region N atds 


PG1528+062 


3 


PG0918+029 


5 


PG1530+057 


3 


PG0942-029 


5 


PG1633+099 


5 


PG1047+003 


4 


PG1657+078 


4 


PG1323-086 


5 


PG2213-006 


4 


PG1528+062 


3 


PG2331+055 


3 


SA 95 


5 


MARK A 


4 


SA104 


4 






SA107 


4 


TOTAL 


26 


TOTAL 


35 



Usually one or more Landolt fields were re- 
observed nightly with an air-mass range of at least 
0.70 in order to measure the coefficients of the atmo- 
spheric extinction of the SPM site, which has excel- 
lent sky conditions. To improve the accuracy, preci- 
sion, and efficiency of the photometric observations 
when required, the filters were observed in forward 
and backward sequences (i.e. UBVRI-IRVBU) , espe- 
cially for the large air-mass observations. 

2.2. Data reduction 

The usual (night and run) calibrations for 
CCD photometry were done during each of our 
observing periods (i.e. bias, twilight-sky flat fields, 
and dark-current determinations) to determine the 
(night and run) mean correcting frames. Standard 
data reduction procedures have been used within 
IRAFE the CCDRED and DA OPHOT tasks (aper- 
ture and PSF photo metry, see iHowell I ll989L Il990t 
Stetson I [l987l Il990h . More details concerning the 
instrumentation and the observing and reduction 
procedures of this project will be given in the 
near future in the su cceeding paper of this project 
( Michel et al. 2010l and references therein). To 
obtain the magnitudes and colors on the standard 
system for the stars associated with these clusters, 
we followed Jordi et al. (1995), and Rossello et 
al. (1998, and references therein). We proceeded 
twofold: 



S IRAF is distributed by NOAO (operated by the Associ- 
ation of Universities for Research in Astronomy, Inc.) under 
cooperative agreement with NSF 



4 



i. AKKAYA ET AL. 



i) The natural magnitude of the filter N is defined as: 
Aat„ = —2.5 • log (ADU's)n, where Xn stands for 
the corresponding filters U,B,V,R, and /, ADU's 
for the analog-to-digital counts, and the subscript n 
for the corresponding quantity in the natural pho- 
tometric system. The atmospheric extinction co- 
efficients for a given filter have been estimated by 
transforming the nightly Xnu's to the corresponding 
magnitudes in the standard system, Ajv s 's, with the 
following equation: 

Xn s — = (zero point) N — k n -X Nii — k n 12 -X Nn -(Ai- 

(1) 

where X^ n is the air-mass when measuring Xmu- 
The subscript N, 12 indicates that the color (Ai — A 2 ) s 
was used to determine the second-order extinction 
coefficient of filter N. Here we follow the conven- 
tion that the effective wavelength \\ e ff < Meff to 
construct the color (Ai — A2). Finally, for a proper 
determination of k n and k n 12 by a least square 
solution, sufficiently large ranges in the airmasses 
and colors of the standard stars (AXjv > 0.7 and 
A(Ai-A 2 )„ > 0.8 for SPM) must be obtained. Note 
that the standard magnitudes and colors are known 
to an accuracy of about two percent, reflected in the 
errors of the final transformations of the (bright) 
standard stars, and that the observed magnitudes 
and airmasses are measured quantities that can have 
an even better precision. To further simplify the 
equations, the extra-atmospheric instrumental mag- 
nitudes were then introduced using the extinction 
coefficients of the night: 

AjVi = Ajv„ — K N -Xjy n — K N12 -Xjy n -(\i— A 2 ) s . (2) 

An instrumental color is the subtraction of two 
instrumental magnitudes with different passbands, 
(Ai — A 2 )i = Xn — A 2 i. 

ii) Once the atmospheric extinction coefficients k n 
and k n 12 have been determined and applied, the 
nightly transformation coefficients are calculated 
(i.e. /?i2 and 712) with the following relations for 
the colors: 

(Al-A 2 )i = K0,12+y<3l2-(Al-A2) s +7l2-(Al-A2)s (3) 

Due to the Balmer discontinuity that lies in both 
the U and B passbands, a better transformation for 
the U-B color has been achieved by substituting the 
quadratic term on the right side of the above equa- 
tion with a linear term in the color B-V, obtaining 
the following expression: 

(Ai - A 2 )j = K ,i2 +(3i2 ■ (Ai - A 2 ) s +712 • (A 2 - X 3 ) s , 

(4) 



where Ai e // < A2 e // < X^eff- For the case of the 
magnitude V, Equation (3) has been used as follows: 

Vi-V s = koi + Pu • (Ai - A 2 ) s + 712 • (Ai - X 2 f s . (5) 

For Equations (3)-(5), ^0,12 and «oi & r e the zero- 
points of the transformations of the colors (Ai — A2) s , 
i.e. U-B, B-V, V-R, V-I, and of the V magni- 
tude, respectively. The coefficients /3i2 and 712 are 
the respective first- and second-order transformation 
coefficients. 

A2)s In general, the second-order atmospheric extinc- 
tion coefficient n VR is expected to be close to zero 
due to the nearly constant level (ozone-band contri- 
bution) of the atmospheric extinction curve at SPM 
near 5500A (Schuster & Parrao 2001). The second- 
order extinction and linear-transformation coeffi- 
cients for correcting to extra-atmospheric standard 
magnitudes and colors are very similar from night 
to night, and also from run to run, because, i) the 
SPM has excellent sky conditions, and ii) the same 
instrumental setup, observing techniques, and data 
reduction procedures were used for all nights during 
both observing runs. In Table 2 the mean zero-point 
corrections, atmospheric extinction, and transforma- 
tion coefficients are given. 

In Tables 3, 4, and 5 are given the final 
transformed CCD UBVRI photometric values 
for the open clusters, Be 89, Ru 135, and 
Be 10, respectively. In the text only the first 
three lines of each table are shown as exam- 
ples. The full versions will be available online 
at the CDS and WEBDA. In these tables the 
columns 1 and 2 give the following: X and 
Y (pixels), the position of a star in the CCD 
field; columns 3, 5, 7, 9, and 11: the magni- 
tude and color indices, V, (B-V), (U-B), (V- 
R), and (V—I), respectively (in magnitudes); 
and columns 4, 6, 8, 10, and 12: the respective 
photometric errors, oy, as-V: &XJ-B, o~v-R> an d 
oy_j (in magnitudes), as provided by IRAF. 

2.3. The data inspection tools ELIPSE and SAFE 

Since the stellar density of a cluster increases to- 
wards its center with respect to the fi eld stars, an 
AWK macro (ELIPSE, iMoitinho lliool was used to 
extract the data of the central region of a given clus- 
ter, as defined by visual inspection in a visual (V) 
or red (R) image, thus increasing the contrast of the 
cluster with respect to the surrounding field stars. 
An ellipse was fitted visually to the image in order 
to extract the photometric data of the central region 
of the cluster. To further support the analyses of 
the clusters, a java-based computer program (SAFE, 
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TABLE 2 

ATMOSPHERIC EXTINCTION AND TRANSFORMATION COEFFICIENTS 



color Ai A 2 A 3 k 1 K i,i2 ^0,12 P12 712 rms 



July 2001 



(U-B) 


UBV 


0.472 


-0.056 


+ 1.625 


0.711 


+0.263 


0.028 


(B-V) 


BV - 


0.243 


-0.050 


+0.409 


1.016 


-0.050 


0.010 


V 


V R- 


0.106 


+0.079 


+2.375 


0.033 


-0.008 


0.016 


(V-R) 


V R- 


0.104* 


+0.030* 


+0.027 


0.973 


+0.011 


0.012 


(V-I) 


VI — 


0.087* 


-0.035* 


-0.151 


0.923 


+0.070 


0.010 








February 


2002 








(U-B) 


UBV 


0.325 


-0.056 


+ 1.765 


0.751 


+0.313 


0.037 


(B-V) 


BV - 


0.212 


-0.050 


+0.470 


0.979 


-0.023 


0.018 


V 


VR- 


0.082 


+0.079 


+2.455 


0.035 


-0.054 


0.027 


(V-R) 


VR- 


0.054* 


+0.030* 


-0.000 


1.023 


-0.008 


0.012 


(V-I) 


VI- 


0.056* 


-0.035* 


-0.165 


1.038 


+0.004 


0.014 



* indicates that extinction coefficients refer to A2, otherwise to Ai 



TABLE 3 



CCD UBVRI photometry of Be 89 

X Y V oy (B-V) a B-v (U-B) Qu-b (V-R.) ay-R (V-I) <Jy~ 

767.7 538.3 11.261 0.006 0.441 0.009 0.002 0.007 0.315 0.009 9.999 9.999 

466.5 496.2 12.118 0.011 1.362 0.018 1.191 0.011 0.742 0.013 1.399 0.013 

511.7 742.6 12.172 0.008 0.443 0.011 0.268 0.005 0.296 0.011 0.569 0.017 



Note to table: The full version of this table is available online at the CDS and WEBDA. 



TABLE 4 



CCD UBVRI photometry of Ru 135. 

X Y V cry (B-V) a B-v (U-B) <tu-b (V-R) cry -a (V-I) Qy -i 

862.7 231.5 11.132 0.003 0.597 0.004 0.212 0.002 0.360 0.010 9.999 9.999 

237.1 325.4 11.689 0.003 0.710 0.003 0.364 0.003 0.387 0.014 0.860 0.011 

294.0 387.3 11.871 0.008 0.742 0.003 0.405 0.004 0.450 0.038 0.895 0.014 



Note to table: The full version of this table is available online at the CDS and WEBDA. 



TABLE 5 



CCD UBVRI photometry of Be 10. 

X Y V ay (B-V) a B-v (U-B) au- B (V-R) cry-R (V-I) ay-i 

374.5 153.4 11.244 0.004 0.433 0.004 0.246 0.002 0.220 0.005 0.500 0.004 

528.6 217.9 12.352 0.002 1.483 0.004 1.346 0.009 0.851 0.003 1.631 0.003 
804.6 881.2 12.746 0.003 0.738 0.014 0.231 0.004 0.440 0.014 0.873 0.005 



Note to table: The full version of this table is available online at the CDS and WEBDA. 



McFarIandll2010( l was developed and used to help 
us in the vis ualization and analysi s of the photomet- 



ric data fe.g. lSchuster et al. M2007T ) . These programs 



facilitate the elimination of field and apparent non- 



member stars of a given cluster from the diagnostic 
diagrams used to enhance the apperception of clus- 
ter features. Once a satisfactory first estimate of the 
parameters was obtained, a full-frame solution was 
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also consulted and refined. 

SAFE is capable of displaying simultaneously in 
different color-color (CC) and color-magnitude (CM) 
diagrams the cluster's data and has an interactive 
way to identify a (group of) star(s) in one partic- 
ular diagram and to see where it falls in the other 
diagrams. This program is capable of displaying up 
to 16 different diagrams for a given cluster and is 
very useful for the determination of a cluster's phys- 
ical parameters. Figure l(a-c) presents the DSS red- 
filter images of Be 89 (panel a), Ru 135 (b) and Be 10 
(c), with the regions analyzed in this work enclosed 
by ellipses. The central (X, Y) pixel coordinates of 
the nearly circular regions in Figure l(a-c), which 
are considered for the photometric analyses are the 
following: (584, 488) pixels for Be 89, (542, 504) for 
Ru 135, and (517, 493) for Be 10. The diameters 
in arcminutes (AX, AY) of nearly circular regions 
in Figure l(a-c) are the following: (2.27, 2.65) for 
Be 89, (2.62, 2.67) for Ru 135, and (3.12, 2.34) for 
Be 10. 

3. ANALYSES OF THE OPEN CLUSTERS 
BE 89, RU 135, AND BE 10 

The (U-B, B-V), two-color or CC, diagram, and 
five CM diagrams have been used together with the 
zero-age-main-sequence (ZAMS) intrinsic-color cal- 
ibrations of Schmidt-Kaler (1982, hereafter SK82) 
and with the Padova isochrones (Girardi et al. 2000, 
hereafter GBBC; Bertelli et al. 2008; Marigo et 
al. 2008, MGBG) to obtain reddenings, metallicities, 
distance moduli, and ages for these clusters. 

Our analysis technique for our program clusters 
places particular emphasis upon the fit of the ZAMS 
intrinsic colors and Padova isochrones to the obser- 
vational data of the clusters, and this depends in turn 
upon important characteristics of the CM and CC 
diagr ams for open clusters (e.g. iPaunzen fc Netopil 
2006l their §3), which are summarized as follows: 



1. A procedure for eliminating non- members, 

2. A determination of the interstellar reddening as 
accurately as possible, 

3. Visibility of the turn-off, 

4. Compensation for binary stars which tend to 
widen the main-sequence distribution, 

5. Consideration of the red-clump stars (if present) 
to improve the isochrone fit, and 

6. An appropriate choice of the isochrone which 
corresponds to the correct heavy-element abun- 
dance (Z). 



Regarding the locus of the main sequence in a 
CM diagram, and independent of any cosmic disper- 
sion, the main-sequence strips or bands in the CM 
diagrams are affected by the contamination of binary 
and multiple stars; particularly, the mid-points are 
shifted to brighter magnitudes and the colors to red- 
der values due to this contamination, and also some- 
times due to variable intercluster extinction. For this 
reason, the SK82 ZAMS and the MGBG isochrones 
have been fitted to the blue- and faint-most concen- 
trations within the observed broad main-sequence 
bands whenever possible, assuming that these con- 
centrations re flect the single-star distribution (e.g. 



Carney 1 12001|) . and that most stars observed red- 



and bright-ward of these are in fact binary, or mul- 
tiple, systems. 

In the absence of proper- motion/radial- velocity 
measurements to insure cluster membership, and to 
minimize the effects of field-star contamination, we 
have concentrated more on the central regions of 
the clusters rather than using the full-frame CCD 
images; this has been accomplished by using the 
ELIPSE or SAFE programs, described above. These 
have been used to select an elliptical, or polygonal 
area (with as many as 10 sides), centered on the open 
cluster as seen in a V or R image, excluding stars 
outside this area from further analyses. (See Fig- 
ure 1). These interactive analyses increase greatly 
the contrast of cluster members with respect to the 
field stars, and thus the scatter in the CM and CC 
diagrams is significantly reduced. 

Also, the observational errors, e.g. C(u~B)i of 
these three clusters have been considered as a cri- 
terion in selecting the more reliable data for further 
analyses. The values of U(jj_b) are almost always 
larger than the ones of a^-v) due to the smaller 
sensitivity of the CCD in the ultraviolet, and the 
errors <j(R_n, ovy_n, and &(v~r) are among the 
smallest. The observational errors, such as 
tT (B—v)i an( l a {B-R)-< have been selected to be less 
than ps CP 10 (and sometimes < 0™05) in some of 
the diagnostic diagrams presented in the analyses 
to follow, such as the (U-B,B-V), (V,B-V), and 
(V,B-R) diagrams. 

Interstellar reddenings of the program clusters 
have been estimated from shifts of the intrinsic-color 
sequences of SK82 in the (U-B,B-V) diagram, until 
the best fit to the data of the clusters: along the U-B 
axis by 0.72E(B-V)+0.05E(B-Vj 2 and along the (B- 
V) axis by E(B-V). For this, F-type stars have been 
fitted above the main sequence of SK82 (i.e. blue- 
ward in (U-B)), and simultaneously the red-clump 
stars above the red-giant colors of SK82 with consis- 



UBVRI PHOTOMETRY OF OPEN CLUSTERS 



7 



tent ultraviolet excesses according to the normaliza- 
tions of Sandage (1969). The two-color sequence of 
SK82 has been constructed from the intrinsic colors 
of SK82 for zero-age dwarfs ((B-V) < m 75) and 
for giants ((-B-V)o > 0'. n 75). 

Once the two-color sequence of SK82 has 
been fixed as indicated above, to determine the 
photometric metal abundance, [Fe/H], one first 
locates the F-type stars in the (U-B, B-V) diagram 
and compares their location with that of their 
counterparts of known metallicity (e.g. the SK82's 
ZAMS calibration). Deviations between the two are 
due mainly to their differences in metal content, 
an ultraviolet excess, 5(U-B), being caused by 
differences in line blanketing. The metal-deficient 
F-type cluster stars, if present, lie blue-ward of the 
"hump region" of the ZAMS sequence, where an 
eyeball-fitted osculating curve similar to "the hump" 
has been fitted to the data points of the F-stars 
(i.e. the thick line above the hump of the F-star 
region in Figures 2, 5, and 8) and, simultaneously, 
to the red-clump stars (if present), since they 
also will lie blue-ward of the red-giant colors of 
SK82 with a corresponding ultraviolet excess. This 
ultraviolet excess is correlated with the photometric 
metallicity of the cluster. Then, a metallicity 
value, [Fe/H], for a cluster can be derived from the 
empirical calibration, [Fe/H]-5({/-B)o.6, of Karatas 
& Schuster (2006), allowing the determination of 
[Fe/H] independently of the isochrones to be fitted 
to the data, thus reducing from three to two the free 
parameters to be derived from the CM diagrams. 
Heavy-element abundances (Z) of the three clusters 
have been obtained from the photometric metal 
abundances [Fe/H] with the expression 



10^1. Z, 



-0.019. 



(0) 



Finally, the appropriate isochrones of MGBG were 
computed online in terms of the resulting heavy- 
element abundance for further analyses of the clus- 
ters (distances and ages). 

To estimate the the age of a cluster (A) and 
the true distance modulus (DM= Vq-Mv) in a 
CM diagram, for example the (V, B-V) diagram, the 
absolute magnitudes, My , of the MGBG isochrones 
have been shifted by DM + 3.LE (B - V) along the 
magnitude axis and their corresponding colors, (B— 
V)q, reddened by adding the color excess E(B-V) 
until some DM value provides a good fit of the ap- 
propriate isochrone to the faint/blue concentration 
of the observed lower main sequence of the clus- 
ter and, if present, of the red-clump stars. One 



has to take into account when determining the DM 
that metal-poor stars are sub-luminous as compared 
with their solar-like counterparts by determining a 
reliable value for Z from the CC diagram. To in- 
fer the age of the cluster, the (logarithm of the) 
age of the isochrones, log A, has been varied until 
a good match with the observed sequences, i.e. the 
upper main-sequence (MS), the turn-off (TO) stars, 
and, if present, the red-clump (RC) stars, has been 
achieved. A fine tuning of the DM has been made 
if necessary. The uncertainties of E(B-V), Z, DM, 
and log A are discussed in Section 4. 

Following a similar procedure to that outlined 
above, the distance moduli and cluster ages have 
also been derived from analyses of four other CM di- 
agrams for each of the clusters. The corresponding 
color excesses applied in the diagrams were iterated 
starting with the previously derived color excess esti- 
mates, E(B-V), and the results were inter-compared 
by means of the standard interstellar e xtinction law 
adopted (see Table 6 below; also cf. iDean et al. I 
ll978HMathi71ll99dlStoizvs lll995l) until satisfactory 
solutions were obtained for all the CM diagrams. 
The derived extinction laws do not differ significantly 
from that of Table 6. 

TABLE 6 

Adopted interstellar extinction law 



E(V-I) 


E(R-I) 


e(v-r) 


E(B-R) 


E(B-V) 


E(B-V) 


E(B-V) 


E(B-V) 


1.25 


0.69 


0.56 


1.56 



3.1. Be 89 

The (U-B, B-V) diagram of Be 89 is shown in 
Figure 2. An interstellar reddening of (B-V) = 
0™60 ± 0™09 has been derived by shifting the in- 
trinsic two-color stellar sequence of SK82 along the 
reddening vector as described in the previous sec- 
tion. (Another possibility, to fit the stars by E(B- 
V) ~ 0™73 to the blue (B-star) branch of the ZAMS 
curve, would leave many stars far from a good fit). 
Six stars apparently in the cluster are noticed with 
(B-V) « 17>6 and (U-B) « 1™4 (big open circles in 
Figure 2) lying near, but above (i.e. blue-ward of) 
the giant sequence of SK82, the expected location of 
the RC stars; their subsequent locations in the CM 
diagrams confirm this classification (see Figures 3 
and 4). A seventh candidate falls further from the 
expected RC locus in four of the five CM diagrams. 

The F-type and RC stars of Be 89 (cf. Figure 2, 
(B-V) rj 1™0 and rs 1™6, respectively) lie above the 
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Fig. 2. The (U-B, B-V) diagram of Be 89. The "S" 
curves (upper parts, ZAMS, and lower parts, red giants) 
have been taken from the two-color relations of SK82 
and are displayed for the interstellar reddening values 
E(B-V) = (TOO and (T60 (the bluer and redder ver- 
sions, respectively). A reddening vector is also shown as 
an arrow, and big open circles mark the six RC candi- 
dates, and open squares, the blue-straggler ones. 
A heavy solid curve represents our best fit to the data 
which has been adjusted to the main-sequence F-type 
stars above (i.e. blue-ward of) the ZAMS colors of SK82 
and, simultaneously, to the RC stars above the red-giant 
colors of SK82. This fit has been used to estimate the 
heavy-element abundance of the cluster, which is shown 
in Table 7. 

(reddened) ZAMS two-color calibration of SK82 by 
S(U-B) ~ (Tl. Our best eyeball fit to the data is 
shown as the heavy solid curve in Figure 2. In the 
dereddened two-color diagram, the heavy line gives a 
value of (U-B)q = -0™10±0 m 02 at (B-V) = (P44, 
and at this same color index, the highest point of 
the SK82 hump has (U-B) = -0™02. The result- 
ing ultraviolet excess, S(U~B) = +0™08 ± (T02, has 
been converted to 6(0.6) = +0 m 10 ± m 02 at (B- 
V) = +0'. n 60 with the normalization ratios given 
by Sandage (1969, his Table 1A). These values have 
been listed in Table 7, together with the correspond- 
ing photometric metallicity [Fc/H]= — 0.35±0.02 dex 
derived with help of the calibration [Fe/H]-<5(0.6) of 
Karatas, & Schuster (2006). Note that the 5(0.6) 
in the notation of the latter authors corresponds to 
delta(0.6) in the notation of Sandage (1969). Apply- 
ing the above relation between [Fe/H] and Z, where 
[Fe/H] has been estimated as —0.35 ±0.02 dex, gives 
Z = +0.008 ± 0.0003. The online isochrones of 
MGBG have been iterated using this metal abun- 
dance when further analyzing Be 89. 



In Figures 3 and 4, the isochrones of MGBG for 
Z = +0.008 have been over-plotted in five CM dia- 
grams: V, (B-V), V, (R-I), V, (V-I), V, (V-R), and 
V, (B-R) after reddening the isochrones along the 
color axis with a color excess corresponding to E(B- 
V) = 0™60, converted with help of Table 6, and 
adding a visual extinction of Ay = 3.1 • E(B-V) = 
1™86 to the absolute magnitudes of the isochrones. 
The isochrones have then been shifted vertically to 



11 
12 
13 



: (a) 



14 r 

15 r 

16 r 

17 1 



18 



19 - 

1 ^ - E(B-V>= 60 




Z(iso) = 0.O08 
log(A)= 9.45-9 55 
V-M.,= 1 1 .9 



: (b) 



20 



11 
12 

13 7 

14 : " ■ «. 

15 - 

16 h 

17 r 

18 : 

Z(iso) = 0.008 ... r . 

log(A)= 9.50-9.60 ;'.VJ'.j' - 

A Q _ V-M v = 1 1 .9 •'' 
' ^ - E(R-I)= 0.41 
" 0.03 





20 



0.3 



0.6 



0.9 



1.2 



(R-i) 



Fig. 3. CM diagrams, (V, B-V) and (V,R-I), for Be 89. 
Solid lines show the interpolated Z = +0.008 isochrones 
of MGBG (cf. Table 7 for the inferred metallicity). Big 
open circles denote the RC candidates, and open squares, 
the blue-straggler ones. See the text and Table 8 for the 
inferred values of the distance modulus and age. 
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Fig. 4. CM diagrams, (V, V-I), (V, V-R), and (V, B-R), 
(top, middle and bottom panels, respectively) for Be 89. 
The isochrone curves and the symbols have the same 
meaning as in Figure 3. See the text and Table 7 for the 
inferred values of reddening and metallicity, and Table 8 
for the distance modulus and age. 
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Fig. 5. The CC diagram of Rul35. The SK82 curves 
and the symbols have the same meaning as in Figure 2. 
See the text and Table 7 for the inferred values of the 
reddening and metallicity. 

obtain the best fit to the observed lower-MS and and 
RC sequences. This vertical shift is the (true) dis- 
tance modulus, DM = (V -M v ). The best fit for 
Be 89 is DM = ll m 90±0 m 06 (d= 2.4 ±0.06 kpc, cf. 
Table 8). 

To derive an age estimate for Be 89, the 
isochrones of MGBG for Z = +0.008 have been 
shifted in the CM planes as above, i.e. Mv+SAE(B- 
V) + DM and C* (Ai - A 2 ) + E[C(Xi ~ A 2 )], respec- 
tively, where the latter color excesses have been com- 
puted from E(B-V) with help of Table 6, and then 
the isochrone ages have been varied until a satis- 
factory fit to the data has been obtained through 
the observed upper-MS, TO, and RC sequences of 
the cluster (cf. Figures 3-4). The resulting inferred 
mean age of Be 89 is log(A) = 9.58 ± 0.06 dex 
(A = 3.8 ±0.6 Gyr). 

For all of these CM diagrams of Be 89, two 
isochrones have been plotted to provide a means 
for appreciating the uncertainties of the derived dis- 
tances and ages. In Table 8 the range in ages pro- 
vided by these isochrone pairs, the final values for 
the distances and ages from each CM diagram, and 
the mean values for each cluster are given. Error 
estimates of (Vb-My) and log(A) are discussed in 
§3.4 below, and the mean results given in Table 8 
have been calculated with Equations (8) and (9) in- 
serting the corresponding parameters summarized in 
the table. 
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3.2. Ru 135 
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Fig. 6. The (V, B-V) and (V, R-I) diagrams for Ru 135. 
Solid lines show the isochrones of MGBG interpolated to 
Z = +0.004. See the text, and Tables 7 and 8, for the 
inferred values of reddening, metallicity, distance modu- 
lus, and age. Stars shown with open-square symbols are 
most likely field, or blue-straggler, stars. 



The same procedures outlined in §3, and §3.1 for 
Be 89, have also been used for the clusters Ru 135 
and Be 10. A reddening of E(B-V) = 0™63 ± (T12 
has been derived for Ru 135 (cf. Figure 5). How- 
ever, a clump of A-type stars at (B-V) ~ 0^8 and 
(U-B) ~ 0™4 seems to be present, with a horizontal- 
like distribution which does not fit satisfactorily the 
reddened two-color ZAMS curve of SK82. These 
stars (Sp A-types) are probably less reddened than 
Ru 135 by ~ 0^3 in E(B-V), nearer, and most prob- 
ably not cluster members (cf. the open squares in the 
CC and CM diagrams of Figures 5, 6, and 7), or they 
could be blue stragglers belonging to the cluster. For 
this latter case, they would be peculiar because of an 
ultraviolet-flux excess present in their spectral en- 
ergy distributions (SEDs), and only a spectroscopic 
study with good signal-to-noise ratios would reveal 
more about their true nature. 

Ru 135 contains a considerable number of F- and 
later-type stars, and appears to have its blue-most 
turn-off limit at (B-V) « +l m and (U-B) « 0™4, 
corresponding to a dereddened (B-V) ~ +0™43 (i.e. 
Sp « F5V) . The best fit to the observed F-hump se- 
quence in the (U-B, B-V) diagram is the solid curve 
shifted blue- ward with respect to the two-color SK82 
curve (cf. Figure 5). From the ultraviolet excess of 
these cluster F stars and following the procedure out- 
lined at the beginning of §3, [Fc/H]= -0.71 ± 0.02 
dex (Z = +0.004 ± 0.0002) has been derived. The 
isochrones of MGBG with this metallicity have been 
computed on line and used in the following analyses. 

The five CM diagrams, (V, B-V) through (V, B- 
R), of Ru 135 are displayed in Figures 6 and 7 to- 
gether with the reddened isochrones that best fit 
the data for the derived color excess and metallic- 
ity, E(B-V) = m 63 and Z = +0.004. The distance 
moduli, (Vq-Mv), and ages, A, found from these five 
CM diagrams and their respective isochrone fittings 
are given in Table 8. 

In these CM diagrams a significant number of 
stars are seen extending to brighter magnitudes and 
red-ward from the fainter and redder observational 
limits of the main sequences, i.e. the stars extending 
red-ward and upward from (V, B-V) « (18'. n 5, l m 5) 
or V(R-I) w (18 m 5,0™9) (cf. Figure 6). These 
are probably field red-giant stars contributed by the 
Galactic bulge, as suggested by the Galactic longi- 
tude and latitude of Ru 135, I c± 16-4° and b ~ 
+6.2 ° (see lBinnev fc Merrifield"lll998t Istanek et al 



1996l Figures 3.5 and 2, respectively). The fact that 
Ru 135 lies near the direction of the Galactic cen- 
tral region also explains the significant number of 
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Fig. 7. The (V, V-I) , (V, V-R) and (V, S-JJ) diagrams 
for Ru 135. The symbols are the same as in Figure 6. 
See the text, and Tables 7 and 8, for the inferred values 
of reddening, metallicity, distance modulus, and age. 

brighter and bluer foreground stars seen in its CC 
and CM diagrams. 

3.3. Be 10 

In Figure 8 the loci of stars observed in the direc- 
tion of Be 10 are shown in the (U-B, B-V) diagram, 
together with the standard interstellar reddening 
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Fig. 8. The (U-B, B-V) plot of Be 10. The symbols and 
the curves are the same as in Figure 2. 

vector and the two-color curve of SK82, shifted along 
this vector to procure the best fit to the data. From 
the fits along the (B-V) and (U-B) axes, E(B- 
V) = (F75 ± m 09 and ([Fe/H], Z) = (-0.49 ± 0.02 
dex, +0.006 ± 0.0003) are found, following the pro- 
cedures described in §3 and §3.1 (see Table 8 for 
partial and mean results). Again, the appropriate 
isochran.es of MGBG have been computed online 
with this corresponding metallicity and are used be- 
low for further analyses of Be 10. 

For Be 10, DM = (V - 3.1 • E(B-V) - M v ) = 
11?16 ± m 06, the distance, d = 1.7 ± 0.05 kpc, 
the metallicity, Z = +0.006 ± 0.0003, log(A) = 
9.06 ± 0.05, and the age, A = 1.08 ± 0.08 Gyr have 
been measured. Our results are listed in Tables 7 
and 8. The resulting (best) isochrone fitting to the 
corresponding Be 10 data in the (V, B-V), (V, R-I), 
(V,V-I), (V, V-R) and V,(B-R) diagrams are dis- 
played in Figures 9 and 10, where one can see that 
the isochrones reproduce well the observed lower and 
upper MS, the TO, and RC sequences of this cluster. 



3.4. Estimated errors and weighted averages 

In Table 7, the ultraviolet excesses and the metal- 
licities are given for Be 89, Ru 135, and Be 10, and 
in Table 8, the mean values for the distance moduli, 
heliocentric distances, logarithmic ages, and ages, to- 
gether with the corresponding estimates of precision. 
The err ors were calculated in a straigh tforward man- 
ner (cf. iBevington fc Robinson" 2003 . and references 



12 



i. AKKAYA ET AL. 



5(a) 



0.5 



11 
12 
13 
14 
15 
16 
17 
18 
19 
20 

11 

: (b) 

12 L ^ ' 

13 r 

14 r 

15 r 

16 r 

17 r 

18 r 

Z(iso) = 0.006 
log(A)= 9 10-9 20 
HQ _ V-M v = 1 1.10 
- E<R-I)= 52 

: "(R-i,- - 01 
20 — — l 




1 1.5 
(B-V) 



2.5 




0.5 



1 

(R-l) 



1.5 



Fig. 9. The (V, B-V) and (V, R-I) diagrams for Be 10. 
Solid lines show the isochrones of MGBG interpolated to 
Z = +0.006. The larger open circles identify the RC can- 
didates. See the text, and Tables 7 and 8, for the inferred 
values for reddening, metallicity, distance modulus, and 
age. 
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Fig. 10. The (V,V-I), (V,V-R) and (V, B-R) diagrams 
for Be 10. The symbols are the same as in Figure 9. See 
the text, and Tables 7 and 8, for the inferred values for 
reddening, metallicity, distance modulus, and age. 

therein). In the following the details of this error 
analysis are presented. 

3.4.1. Errors in E(B-V) and Z 

The random errors in the color excess E(B-V) 
and photometric metallicity [Fe/H] were estimated 



UBVRI PHOTOMETRY OF OPEN CLUSTERS 



13 



follows: 

) By moving the two-color curve of SK82 back- 
ward and forward along the standard reddening 
vector in the (U-B, B-V) diagram until a good 
fit with the observed GK-type, RC, and F-hump 
sequences was achieved. (No BA-type sequences 
were present for these clusters.) The precision 
of the determinations also depends on the scat- 
ter of the data points (cf. Figure 5 of Ru 135 and 
Figure 8 of Be 10 for good and regular cases, re- 
spectively). The uncertainties given in Tables 7 
and 8 reflect this. Following this procedure, a 
typical error in E(B-V) for the quality of the 
reduced data of our clusters is (conservatively) 
« 07 1 04. The systematic error in E(B-V) de- 
pends on the color calibration used. In the case 
of SK82, the uncertainty can be safely assumed 
to be, at the most, of the order of the difference 
between two adjacent spectral subclasses. 

) The random photometric-metallicity uncer- 
tainty has then been estimated from the 
parabolic (eyeball) fit to the data of the maxi- 
mum characterizing the ultraviolet flux excess 
of the stars at the dereddened color (B-V)o 
~ 0^44 (Sp ~ F5) and then following Sandage's 
(1969) normalization procedure. 

The uncertainty of the metal content Z 
was determined from th e relation (e.g. 
Bevington fc Robinson | 2003): 



<Tz = In 10 X Z X <T[Fe/H]- 



(7) 



<J[Fe/H] nas been estimated from the uncer- 
tainty in the ultraviolet excess S(U-B) at the F 
hump between the observed and the SK82 two- 
color curves and is typically ±0™02. Assuming 
(Z) = 0.006 (the mean of the three clusters) 
| az |< 0.0003 is obtained with Equation (7) 
above. Assuming an error of about 0.001 for 
Z is, in our case, a quite conservative estimate. 



) On the other hand, the deviation of the as- 
sumed reddening vector from the "true" one 
depends on the quotient , which can 

strongly deviate l ocally from its canonical 



value of 0.72 (see IChavarria-K et al. I 11987 



Johnson Ill977l and references therein). This un- 



certainty may produce errors larger than the 
precision quoted above. For a crude estimate, 
using the extremes of the cited values of ^fB-v) 
and a typical color excess of E(B-V) = 0™50, 



the uncertainty in 6(U-B) could be as large as 
0™15-0™20. However, since our displacements 
of the SK82 curve in the CM-diagrams are con- 
sistent with the canonical value for the inter- 
stellar extinction law, we assume that this error 
contribution is negligible statistically. 

iv) Another systematic uncertainty results from the 
two-color calibration of SK82: from the uncer- 
tainty of (U-B) or (B-V) in the two-color cal- 
ibration curve of SK82, which is expected to 
be of the order of the difference between two 
typical spectral subclasses (in our case < 0™05) 
for the (U-B) index; and from the fit of the 
whole curve to a cluster data set, of the order 



m 05 



where there are N pivotal points con- 
sidered when adjusting the SK2 two-color curve 
(the BA-type, the F-hump, the GK-type, and 
the RC sequences; i.e. N > 3). Summarizing, 
the systematic uncertainty in S(U-B) should be 
less than roughly three times that given by the 
precision of the flux measurements. 



3.4.2. Errors in the distance moduli and ages 

i) The uncertainties o~dm% for the moduli given in 
Table 8 that result from fitting the appropri- 
ate isochrones to the data in the CM diagrams 
depend on the photometric uncertainties (flux 
measurement and standard-transformation er- 
rors) , the absolute magnitude and intrinsic color 
calibration errors (see for example, SK82), the 
color excess uncertainty of a given color (which 
depends on E(B—V)), and on the reddening law 
adopted. We assume that the isochrones only 
contain the errors of the absolute magnitude and 
intrinsic-color calibrations and that the pho- 
tometric and transformation errors are small 
(« 0™03) when compared to the other sources 
of error. In our case, the largest contribution to 
the distance modulus uncertainty is due to the 
uncertainty in the absolute-magnitude scale, fol- 
lowed by the uncertainty in the slope of the red- 
dening vector, and the color-excess error, about 
0™3, 0™15 and 0™12, respectively, which com- 
bine to give an expected total uncertainty as 
large as an ah — 0". n 25. 

The moduli resulting from the CM diagrams of 
each object and the mean moduli for the three 
clusters are given in Table 8, and the mean of 
the moduli has been derived from the five mod- 
uli, weighted with their respective (usually un- 
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equal) precisions, with the following expression: 
ThMi -; (* = !,..., 5), (8) 



DM 



Z(DM) l /a 2 L 



and the associated mean uncertainty is esti- 
mated from the individual uncertainties of the 
five CM-diagrams of a given cluster by the rela- 
tion: 

The combined error is the square root of the sum 
of the squared uncertainties and is expected to 
be about 0™15, or even less. 

ii) The uncertainty in the log (A) has a random er- 
ror due to the (eyeball) fit of the isochrone with 
the appropriate metallicity to a given CM di- 
agram of a cluster in question, and a quanti- 
tative estimate is obtained by jiggling bright- 
ward and faint-ward the isochrone curve until 
a good fit of the lower main sequence produces 
the DM. Then the age of the isochrone is var- 
ied until a good fit to the upper main sequence, 
the TO, and the RC sequences is achieved. The 
two isochrones shown in the CM diagrams of the 
program clusters give a quantitative estimate of 
this last error. Several different authors have 
computed isochrones as function of the metallic- 
ity, and the physics behind seems to be well un- 
derstood. One does not expect a large variation 
in the log (A) error due to any uncertainty in the 
physics, and the uncertainties of E(B-V) and 
(Vo-My) play a secondary role because the age 
errors depend more on the form of the isochrone 
curve and how it embraces the data (i.e. the the 
upper main-sequence and TO regions and the 
RC sequence) and, less significantly, on the red- 
dening law (except perhaps the blue and near- 
ultraviolet filters) . More problematic is the case 
when the TO region is not well defined (i.e. iso- 
lated from field stars) and/or the RC sequence 
is not present. In our case, the errors for the 
different colors of Table 8 reflect these uncer- 
tainties. 

4. COMPARISON OF FUNDAMENTAL 
PARAMETERS OF THE THREE CLUSTERS 
The reddening values of our three clusters have 
been compared to ones derived from the dust maps 
of Schlegel, Finkbeiner & Davis (1998; hereafter, 
SFD); these are based on the COBE/DIRBE and 
IRAS/ISSA maps, and take into account the dust ab- 
sorption all the way to infinity. E(B-V)(£,b) co val- 
ues of our three clusters have been taken from SFD 



maps using the web pages of NELq These E(B- 
V)(£,b) QO values are m 99 for Be 89, and l m 06 for 
both Ru 135 and Be 10. However, Arce & Good- 
man (1999) caution that SFD reddening maps over- 
estimate the reddening values when the color ex- 
cess E(B-V) is more than w 0™15. For the revi- 
sion of SFD reddening estimates, the equations of 
Bonifacio, Monai & Beers (2000) and Schuster et 
al. (2004) have been adopted. Then the final redden- 
ing, E(B-V)a, for a given star is reduced compared 
to the total reddening E(B-V)(£,b) OCl by a factor 
{1 — exp[— dsmb/H]}, where b, d, and H are the 
Galactic latitude (Column 2 of Table 9), the dis- 
tance from the observer to the object (Column 7 
of Table 9), and the scale height of the dust layer 
in the Galaxy, respectively; here we have assumed 
H = 125 pc (Bonifacio, Monai & Beers 2000). Note 
that Galactic latitudes of our three clusters are less 
than 10°. These reduced final reddenings are E(B- 
V) A = m 54 for Be 89, m 36 for Ru 135, and m 64 
for Be 10. 

For Be 89, our reddening value of 0^60 is in good 
agreement with the value of 0^54 obtained from the 
dust maps of SFD. For Be 10 our reddening value of 
E(B-V) = (F75 is within about la of the value (T64 
derived from the SFD dust maps, and for Ru 135, our 
reddening value of 0™63 differs by about 2a from 
the value of ff. n 36 obtained from these SFD maps. 
These reddening values derived by different methods 
are in reasonable agreement with each other, giving 
confidence to our results. 

As can be seen from the summarized results given 
in Table 9, the reddening value 0™60 found here 
for Be 89 is smaller than the E(B-V) = 1™03 of 
Tadross (2008; hereafter T08a), and than t he E(B- 
V) = l m 05 of (jSubramaniam et al. I I20T0L hereafter 
S10). Our derived distance modulus and distance 
for Be 89, {{V -M v ), d(kpc)) = (ll m 90±0 m 06, 2.4 ± 
0.06), are smaller than the values of (12 m 39,3.00) 
of T08a and larger than the (ll m 54,2.04) of S10. 
Our inferred age (log(A), A(Gyr)) = (9.58,3.8 Gyr) 
for this cluster is considerably older than (8.93, 0.85 
Gyr) given by T08a and larger than the estimate 
(9.02, 1.06 Gyr) by S10. For the analysis of Be 89, 
T08a used JHK photometry and the isochrones of 
Bonatto et al. (2004) with a solar metallicity. This 
is, partially, the origin of the disagreement between 
the two age estimates, since our lower metallicity 
for Be 89 will necessarily lead to a larger age for a 
given TO. Also, most probably, the differences are 
partially due to the different procedural approaches 
for estimating the fundamental parameters; we de- 



3 http: / /nedwww. ipac.caltech.edu/forms/calculator. html 
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rive in a straightforward manner the estimates for 
the interstellar extinction and metallicity: by fitting 
SK82's ZAMS to the data in the (U-B, B-V) di- 
agram, by then measuring the ultraviolet excess of 
the F-type stars to derive a cluster metallicity, and 
finally using the appropriate isochrones in CM dia- 
grams to estimate the true distance modulus and age 
of Be 89. Two parameters (reddening and metallic- 
ity) are estimated in a CC diagram separately from 
the other two parameters (distance and age) from 
the CM diagrams. S10 have also assumed a solar 
metallicity (Z@) for their isochrones (from GBBC) 
and have used only CM diagrams to estimate the 
reddening, distance, and age of Be 89. 

Previous results in the literature for Ru 135 are 
found in the work by Tadross (2008; hereafter T08b), 
and for Be 10 in the papers by Lata et al. (2004; L04) 
and Maciejewski & Niedzielski (2007; MN07). Our 
reddening value E(B-V) = 0™63 ± (T12 for Ru 135 
is significantly smaller compared to the reddening 
value of 1™10 given by T08b. Also, our derived dis- 
tance modulus and distance, ((Vo~M v ), d(kpc)) = 
(9™58, 0.81), for Ru 135 are significantly smaller than 
(11^33,1.85), and our inferred age (log(A),A) = 
(9.58, 3.80 Gyr) is considerably older than (8.70, 0.50 
Gyr), values by T08b. 

In defense of the present results, our value for 
E(B-V) (0™63) falls between the value derived from 
SFD (0 m 36) and the value of T08b (1™10); our value 
is in much better agreement with SFD. T08b used 
the solar-metallicity isochrones of Bonatto, Bica, & 
Girardi (2004), and his results are based on the com- 
parison of isochrones to observed data in the J, (J- 
H) and K 7 (J-K) planes of infrared photometry. 
These differences between our values and those of 
T08b are probably due mainly to the largely differ- 
ent values for the interstellar reddening, but also to 
the difference in the assumed metallicities, to the 
use of different stellar models and isochrone sets, 
which make use of differing input physics and colour- 
temperature transformations, and to distinct photo- 
metric data sets. 

For the Be 10 open cluster, our reddening value 
E(B-V) = 0™75 is in reasonable agreement with 
the value of E(B-V) = 0™87 given by L04, and in 
good agreement with E(B-V) = m 71 by MN07. 
For the metallicity of the Be 10 cluster, L04 adopt 
the Z = +0.008 isochrones of Girardi et al. (2002), 
and MN07 adopt the solar isochrones of Bcrtelli et 
al. (1994). From our two-color diagram, Z — +0.006 
has been derived (see §3.3), which is in agreement, 
within the error bars, with the value of L04. Our 
distance modulus and distance for Be 10, ((Vq- 



M y ),d(kpc))= (11^16,1.70) differ from the values 
(11™8, 2.3) of L04, but very little from the values of 
MN07, (ll m 26, 1.79). Our inferred age (log(A), A) 
= (9.06, 1.08 Gyr) for this cluster disagrees by al- 
most a factor of two (0.5 Gyr) with L04, but is in 
good agreement with MN07, (9.00, 1.00 Gyr). Again, 
our interstellar reddening for Be 10, E(B-V) = 0™75 
falls between the value derived from SFD (0™64) and 
the value 0™87 by L04. 

The age values in Table 9 have been compared 
to ages estimated with the (age, AV) calibration 
given by Carraro & Chiosi (1994; their equation (3)). 
Note that this last calibration does not consider the 
metal abundance of the cluster. Here, AV means 
the magnitude difference between the RC and TO, 
which is well known ge indicator. Both open 

clusters Be 89 and Be 10 have RC candidates (see 
the CM plots for these two clusters, Figures 3-4, 
and Figures 9-10, respectively). TO values occur 
at V « 16 m 5 for Be 89 and V « 14 m 8 for Be 10, 
whereas the RCs occur at V « 15 m 3 and V « 14 m 7, 
respectively. From this age-AU calibration of Car- 
raro & Chiosi (1994), ages have been estimated as 
log(^) = 9.1 (1.3 Gyr) for Be 89 and log (A) = 8.6 
(0.4 Gyr) for the Be 10. 

The average age values given by us, log (A) = 9.58 
(3.8 Gyr) for Be 89 and log(A) = 9.06 (1.08 Gyr) 
for Be 10 are somewhat older than the ones esti- 
mated from this relation of Carraro & Chiosi (1994). 
However, these age differences are at least partially 
explained by the sub-solar metallicities of these two 
clusters ([Fe/H]= -0.35 dcx for Be 89 and -0.49 
dcx for Be 10; see §3.1, §3.3, and Table 9). Lower 
metallicities require larger ages for the same TO. 

In Tabic 9 our results are summarized for Be 89, 
Ru 135, and Be 10: Columns 1 and 2 contain the 
cluster name and Galactic coordinates, respectively. 
The resulting reddening, E(B-V), is given in Col- 
umn 3. The metallicity and heavy-element abun- 
dances, [Fe/H] and (Z), are given in Columns 4 and 
5, respectively. True distance modulus values, (Vo- 
Mv), and their corresponding heliocentric distances 
to the observer are given in Columns 6 and 7, respec- 
tively. Column 8 gives the average age (i.e. log(A); 
where A is in years), as derived from the five CM dia- 
grams. Different isochrones used by us and by other 
authors are referenced in Column 9. Average Galac- 
tocentric distances are listed in Column 10. The cor- 
responding references from the literature are listed in 
Column 11. 

5. CONCLUSIONS 
CCD UBVRI photometry of three poorly stud- 
ied Galactic open clusters, Be 89, Ru 135, and Be 10, 
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has been analyzed, based on new SPM observations. 
The fundamental parameters of reddening, metallic- 
ity, age, and distance of these open clusters have 
been inferred and presented in Tables 7-9. 

The interstellar reddenings and metallicities of 
these three clusters have been determined from two- 
color, (U-B, B-V), diagrams prior to the use of the 
CM diagrams. Heavy element abundances, Z, of the 
three clusters have been found from the ultraviolet 
excess, 5(U—B), of the F-stars by comparison with 
the two-color curve of SK82 (Z$K82 — Zq), by using 
the normalizations of Sandagc (1969), and by ap- 
plying the calibration, [Fe/H]-<5(0.6), of Karatas, & 
Schuster (2006), with the advantage of reducing by 
two the number of free parameters of the isochrones 
when fitting to the data in the CM diagrams. When 
necessary, we have iterated slightly afterwards for 
a better, more consistent, solution for the four clus- 
ter parameters (reddening, metallicity, distance, and 
age). Deeper U frames would improve our deter- 
minations employing this method, which allows us 
to estimate the reddening and metallicity indepen- 
dently using a CC diagram, in contrast to the exclu- 
sive fitting of isochrones to CM diagrams and the use 
of the solar metallicity, which are the more common 
techniques used in the literature. 

The present adjustments of the SK82, CC rela- 
tions to the MS and RC stars, and of the MGBG 
isochrones to MS, TO, and RC stars in the CM di- 
agrams show good consistency and appropriate fits 
for all three open clusters, in the one CC diagram 
and all five CM diagrams. Good consistency is seen 
in the Figures 2-4 for Be 89, Figures 5-7 for Ru 135, 
and Figures 8-10 for Be 10. 

The CC and CM diagrams of Be 89 and Ru 135 
suggest that they are metal-poor and old for their 
location in the Galaxy, compared to other open clus- 
ters. 

For Be 89, stars with V < 16 m 2 and (B-V) < 
0™9 are most likely foreground or blue-straggler 
stars. The blue-straggler and RC candidates in the 
field of Be 89 need spectroscopic and/or astrometric 
observations to test their cluster membership and to 
elucidate their nature. 

Similar candidates for blue-straggler or bright 
foreground stars are seen in the CC and CM dia- 
grams of Ru 135 and Be 10, Figures 5-7 and 8-10, 
respectively. In the case of Ru 135 and for stars with 
V fainter than about 14™2, the onset of the cluster 
in the CM diagrams is clearly seen. Objects brighter 
than this limit and with (B-V) < 0™9 are probably 
blue stragglers or bright foreground stars. 

Despite its similar age to Be 89, no RC stars 



are noticeable in the CM diagrams of Ru 135. On 
the other hand, the CC and CM diagrams of Be 10 
show clear evidence for a RC grouping, although it is 
somewhat younger than the other two clusters. The 
lack of any RC stars in the CM diagrams of Ru 135, 
contrasting with Be 89 and Be 10, may result either 
from relative differences in mass segregation and our 
emphasis on the inner regions of these clusters, or 
from the poorness of these cluster fields and small- 
number statistics. Ru 135, being closer to the Galac- 
tic center, may be more perturbed and less dynami- 
cally relaxed than the other two clusters. Also, Be 89 
and Be 10 each show only eight, or fewer, RC candi- 
date stars, and it is not clear that all of these are in 
fact cluster members. 

For the typical accuracy of photometric observa- 
tions (and we are no exception), the final error esti- 
mates are fixed by the accuracy of the cluster param- 
eters as given by the systematic uncertainties in the 
absolute-magnitude, intrinsic-color, and reddening- 
vector calibrations, for example, the adequacies, or 
not, of the SK82 colors, the MGBG isochrones, and 
the standard interstellar-reddening curve. 

Finally, further radial velocity and proper mo- 
tion information for these clusters will allow us to 
clean with more assurance most non-members from 
the CC and CM diagrams in order to obtain better 
determinations of their physical parameters and to 
better understand the nature of the blue-straggler 
and red-clump candidates in these three open clus- 
ters. Deeper photometric observations, especially in 
the U and B bands, will provide clearer, cleaner, and 
more precise solutions from the CC diagram. 
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TABLE 7 

NORMALIZED (U-B) EXCESSES AND DERIVED METALLICITIES 



Cluster 


(U-B) S K82 


(U-B) ju 


5{U-B) 


J(0.6) 


[Fe/H] 


Z 




[mag] 


[mag] 


[mag] 


[mag] 


[dcx] 




Be 89 


-0.02 


-0.10 


0.08 


0.10 


-0.35 


0.008 


Ru 135 


-0.02 


-0.16 


0.14 


0.14 


-0.71 


0.004 


Be 10 


-0.02 


-0.09 


0.07 


0.11 


-0.49 


0.006 


error 


< ±0.01 


±0.02 


±0.02 


±0.02 


±0.02 


< ±0.001 



TABLE 8 

DISTANCE AND AGE ESTIMATES OF THE CLUSTERS 



Color 


(Vo-Mv) 


d 


log(A) 




A 
SI 




[mag] 


[kpc] 


range 








Be 89: E(B-V)= 0.60 ± 0.09 & Z = 


+0.008 ±0.001 




(B-V) 


11.90 ±0.10 


2.4 ±0.1 


9.45-9.55 


9.55 ±0.15 


3.6+1.4 


(R-I) 


11.90 ±0.20 


2.4 ±0.2 


9.50-9.60 


9.60 ±0.20 


4.0 + 2.3 


(V-I) 


11.90 ±0.15 


2.4 ±0.2 


9.50-9.60 


9.60 ±0.20 


4.0 + 2.3 


(V-R) 


11.90 ±0.10 


2.4 ±0.1 


9.50-9.60 


9.60 ±0.15 


4.0+1.6 


(B-R) 


11.90 ±0.20 


2.4 ±0.2 


9.45-9.55 


9.55 ±0.10 


3.6 + 0.9 


Mean 


11.90 ±0.06 


2.4 ±0.06 




9.58 ±0.06 


3.8 + 0.6 




Ru 135: E(B-V)= 0.63 ±0.12 & Z = 


= +0.004 ±0.001 




(B-V) 


9.50 ±0.15 


0.75 ±0.05 


9.60-9.70 


9.60 ±0.15 


4.0+1.6 


(R-I) 


9.70 ±0.15 


0.87 ±0.06 


9.55-9.65 


9.55 + 0.15 


3.6+1.5 


(V-I) 


9.60 ±0.20 


0.83 ±0.08 


9.55-9.65 


9.55 ±0.15 


3.6 + 1.5 


(V-R) 


9.60 ±0.15 


0.83 ±0.06 


9.55-9.65 


9.60 ±0.15 


4.0 + 1.5 


(B-R) 


9.50 ±0.20 


0.75 ±0.07 


9.60-9.70 


9.60 ±0.15 


4.0 + 1.6 


Mean 


9.58 ±0.07 


0.81 ±0.03 




9.58 ±0.06 


3.8 + 0.7 




Be 10: E(B-A 


/)= 0.75 ± 0.09 & Z = 


+0.006 + 0.001 




(B-V) 


11.20 ±0.11 


1.7±0.1 


9.05-9.15 


9.05 + 0.10 


1.1 + 0.3 


(R-I) 


11.10 ±0.20 


1.7±0.2 


9.10-9.20 


9.10 + 0.10 


1.3 + 0.3 


(V-I) 


11.15 ±0.15 


1.7±0.1 


9.05-9.15 


9.05 + 0.15 


1.1 + 0.3 


(V-R) 


11.15 ±0.20 


1.7±0.2 


9.05-9.15 


9.05 + 0.10 


1.1 + 0.3 


(B-R) 


11.20 ±0.10 


1.7±0.1 


9.05-9.15 


9.05 + 0.05 


1.1 + 0.1 



Mean 11.16 + 0.06 1.70 + 0.05 



9.06 + 0.05 



1.08 + 0.08 
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TABLE 9 



FUNDAMENTAL PARAMETERS OF BE 89, RU 135, AND BE 10 



Cluster 


(Z°, 6°) 


E(B-V) 


[Fc/H] 


Z 


(Vo-Mv) 


d 


log(A) 


Isochrone 


Rgc 


Reference 






[mag] 


[dex] 




[mag] 


[kpc] 






[kpc] 




Be 89 


83.16, +4.82 


0.60 


-0.35 


+0.008 


11.90 


2.40 


9.58 


m8 T 


8.55 


this work 






1.03 




solar 


12.40 


3.00 


8.93 


b4 




Tadross 2008a 






1.05 




solar 


11.54 


2.04 


9.02 


go 




Subramaniam et al. 2010 


Ru 135 


16.42, +6.23 


0.63 


-0.71 


+0.004 


9.58 


0.81 


9.58 


m8 


7.72 


this work 






1.10 




solar 


11.33 


1.85 


8.70 


b4 




Tadross 2008b 


Be 10 


138.62, +8.88 


0.75 


-0.49 


+0.006 


11.16 


1.70 


9.06 


m8 


9.84 


this work 




J 


0.87 




+0.008 


11.80 


2.30 


8.80 


g2 




Lata et al. 2004 




) 


0.71 




solar 


11.26 


1.79 


9.00 


B4 




MN07 



f Isochrone sources: B4 = Bcrtclli et al. (1994); b4 = Bonatto et al. (2004); gO = GBBC; g2 = Girardi et al. (2002); m8 = MGBG 
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Supplementary material 

The following material is available online at the 
CDS and WEBDA: (a) This manuscript, (b) Ta- 
ble 3, which presents standard UBVRI CCD 
photometry and observing errors for the open 
cluster Be 89; the columns 1 and 2 give the 
following: X and Y (pixels), the position of 
a star in the CCD field; columns 3, 5, 7, 

9, and 11: the magnitude and color indices, 
V, (B-V), (U-B), (V-R), and (V-I), respec- 
tively (in magnitudes); and columns 4, 6, 8, 

10, and 12: the respective photometric er- 
rors, cry, <tb-v, &u-b, ov -R-, and oy_j (in mag- 
nitudes), as provided by IRAF. (c) Tables 4 
and 5, which have the same format as Table 3, 
but for Ru 135 and Be 10, respectively. 



